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Abstract 

We investigate quantum gravity in four dimensions using the Regge approach 
on triangulations of the four-torus with general, non-regular incidence matrices. 
We find that the simplicial lattice tends to develop spikes for vertices with low 
^^ coordination numbers even for vanishing gravitational coupling. Dilferent to 

the regular, hypercubic lattices almost exclusively used in previous studies, we 
find now that the observables depend on the measure. Computations with 
nonvanishing gravitational coupling still reveal the existence of a region with 
^-\ well-defined expectation values. However, the phase structure depends on the 

^ii^ triangulation. Even with additional higher-order terms in the action the critical 

[^ behavior of the system changes with varying (local) coordination numbers. 



PACS number(s): 04.60.-Fn, 12.25. 



1 Introduction 

The Regge calculus is an elegant method for investigations of quantum gravity 
through simplicial approximations [1, 2]. The sum over smooth geometries g in 
the Euclidean path integral 

Z= /i?ge-^^(s) (1) 

is replaced by a sum over simplicial lattices with the associated Regge action dis- 
cretizing the gravitational action Ie [3]. Simplicial lattices are specified by their in- 
cidence matrices, describing how the simplices are glued together, and the (squared) 
lengths of the contained links. 

Two complementary approaches have been proposed to evaluate (numerically) 
simplicial path integrals varying either the incidence matrices for equal link lengths 
or, what is called Regge calculus, the link lengths for fixed triangulation. While 
the former method, known as dynamical triangulation, is most successful in two 
dimensions it might have serious problems in four dimensions [4]. The advantage of 
the Regge approach lies in the fact that the classical continuum limit is well-defined 
and quantization in the weak-field approximation can be studied [5]. 

Non-perturbative investigations of the Regge approach in four dimensions have 
been carried out on regular hypercubic triangulations of the 4-torus with rather 
astonishing results [6, 7, 8]: 

(i) Despite the unbounded pure Regge action, a well-defined phase with finite 
expectation values exists for small gravitational coupling parameters [6, 7, 8]. 

(ii) Investigating a family of local measures, their infiuence on the well-defined 
phase was found to be small [8]. 

(iii) Including an additional higher-order B? term in the action, a phase transition 
of 1 order was found, indicating the existence of a non-trivial continuum limit [7]. 

(iv) The incorporation of additional matter fields does not affect the existence 
of the well-defined phase within the Regge approach [9]. 

Until now, the Regge calculus was applied almost exclusively for the regular hy- 
percubic triangulation of the 4-torus. However, if one considers arbitrary topologies 
one has to use necessarily different triangulations and in general one would prefer to 
perform computations on random triangulations [6, 7]. If Regge gravity is a realis- 



tic theory the 'microscopic details' of the underlying triangulation should not affect 
physical expectation values. 

Unfortunately, first studies of the Regge approach for non-regular incidence ma- 
trices signalled possible problems due to the occurence of spikes at vertices with low 
coordination number [10]. To clarify these findings, we investigate in this article 
simplicial quantum gravity on fixed but non-regular triangulations of the 4-torus, 
with the (squared) link lengths as dynamical variables. For this purpose we apply 
barycentric subdivision steps on a regularly triangulated 4-torus, inserting addi- 
tional vertices with small local coordination number. We study the infiuence of 
the measure, the existence of the well-defined phase and the effect of additional R^ 
terms. 

2 Simplicial quantum gravity 

To determine the triangulation of a given 4-manifold Al one has to specify a set 
of N(i d-simplices (i.e. vertices, links, triangles, tetrahedra and 4-simplices) which 
make up a simplicial manifold with the desired topology [3]. Until now, almost 
all numerical computations have been carried out on regular triangulations of the 
4-torus, based on a hypercubic tesselation, with Nq > S"*, Ni = 15iVo, N2 = 50iVo, 
N3 = 60iVo, and N4 = 24iVo. Details of this construction and quantization in the 
weak-field approximation are discussed in the literature [5]. We call a triangulation 
regular if the local coordination numbers are the same for each vertex. In contrast 
to the hypercubic triangulation, we use in this study the construction of Kiihnel 
and Lassmann, triangulating the 4-torus with Nq = 31 vertices [11]. Most interest- 
ingly, the (global) coordination numbers (i.e. Ni/Nq etc.) are the same as for the 
hypercubic case. 

To construct a non-reg(M/ar triangulation we have to change the local coordination 
numbers. We apply a number of barycentric subdivision (b.s.d.) steps, inserting 
additional vertices such that the total number of simplices increases as ANq = 1, 
ANi = 5, AN2 = 10, AN3 = 10, and AN4 = 4. While the number of links attached 
to a vertex is 30 for the regular triangulations it is 5 for the additional b.s.d. vertices 
[10]. 



Once the incidence matrix is fixed one can assign to each link / the squared 
length qi to approximate a smooth 4-geometry. The infinitesimal line element ds^ is 
thus replaced by the finite squared distances {qi} between neighboring vertices. We 
consider the Euclidean sector only, qi > 0, and require that every 4-simplex can be 
embedded in R'^. 

Given a Euclidean configuration {qi} several geometric quantities can be calcu- 
lated. In the following we need the 4- volume Vg of every 4-simplex s and its fatness 
(^s defined as 

(I)s = C 



max/gs(g/4) 



where we set the constant to C = 24. (Cheeger et al. use a slightly different 
but equivalent definition [12].) The fatness takes its maximum value of 0.3125 for 
equilateral simplices and tends to for collapsing ones. Furthermore, we use the 
area At of every triangle t and the associated deficit angle St given by 



St = 27T-J2 






with the dihedral angle O^^t [2, 3, 7]. An important observable is the average curva- 
ture [7] 

The Einstein- Hilbert action is replaced by a sum over all triangles in the lattice 

I d\^R ^ 2Y,AtSt. (5) 

Increasing the number of simplices and decreasing the link lengths such that the 
lattice skeleton approaches a smooth geometry, general relativity is recovered in the 
classical continuum limit if <?!>s > / =const> [12, 13]. 

Concerning quantization, we assign concrete meaning to the functional integral 
(1) replacing the integration over closed 4-geometries by a sum over distinct simpli- 
cial lattices with fixed triangulation. To keep the 4- volume of the simplicial lattice 
finite one incorporates a cosmological- const ant term in the action 

-lE{{qi}) = (3Y,AtSt-^Y.^s. (6) 

t s 

The variables g;. At, Vg are expressed in units of the bare Planck length, the param- 
eter [i determines the gravitational coupling and A fixes the average lattice volume. 



To complete the definition of the simplicial path integral one has to specify the 
integration measure. Within the Regge approach it seems natural to translate the 
gravitational measure into a product over the links of the simplicial lattice 

Dq = Y[fi{l)dqiT{qi,...,qN,), (7) 

where the function T is equal to 1 if the generalized triangle inequalities are fulfilled 
and otherwise, ensuring that the computations are performed in the Euclidean 
sector [3, 6, 7]. The weight function ^u is a (perhaps complicated and non-local) 
function of {qi} associated with each link / [14, 15]. According to previous studies 
we consider as a local and simple ansatz the one-parameter family 



Kl) = qr\ (8) 

where ci > determines the behavior of the measure under rescaling [8]. The case 
(7 = 1 corresponds to the uniform measure and for ci = a scale-invariant measure 
results. 

Considering the behavior of the path integral 



under rescaling one obtains 



^(/3, A) = / n dqiqr'HUl}^ ^' ^'''"' ^^" ""' 



(9) 



-/3{J2At6t) + 2X{J2Vs) = aN^ (10) 

t s 

in the well-defined phase [7]. We set the parameter A equal to a so that (Vg) takes 
the same value Ni/{2N4) at /3 = for all cr > 0. 

Numerical computations within the Regge approach are based on the Metropolis 
method generating a Markov chain of configurations {^/}t according to the measure 
(8) and the action (6) [6, 7, 8]. After thermalization of the system, averages over the 
generated configurations approximate the expectation values. A priori one expects 
that simulations of simplicial quantum gravity with unbounded action might not 
find an equilibrium. However, several studies on regular hypercubic triangulations 
of the 4-torus show a well-defined phase with finite expectation values for (3 < (3c- 
This phase was found to be stable for different measures as long as 5 cr ^ 1 [7, 8]. 



3 Phase structure for non-regular triangulations 

Previous studies indicate that vertices with low coordination number tend to de- 
velop spikes fluctuating with rather long auto-correlation lengths even for vanishing 
gravitational coupling [10]. To analyze the mechanism in more detail, we use a tri- 
angulation as small as possible allowing to perform IM sweeps for each data point. 
We employ the construction of Kiihnel and Lassmann [11] and apply three b.s.d. 
steps to neighboring simplices increasing the total number of vertices to Nq = 34. 
We further use a lower limit for the fatness, <?!>s > / > 0, to restrict the configuration 
space. 

Our first goal is to investigate the infiuence of the measure on the length of the 
spikes for vanishing gravitational coupling, (3 = 0. We performed computations for 
three different values of the measure parameter, a = 1,2'ib' ^^^ decreased the 
limit of the fatness stepwise by a factor i as / = 2™10~^, m = 9, . . . , 0. To study 
explicitly the behavior of the inserted irregularities, we observed the expectation 
value (qi) averaging the (squared) lengths of those links attached to the additionally 
inserted vertices. It measures directly the next-neighbor distances and thus the 
lengths of the developing spikes. This quantity can be compared with the expectation 
value {qj.) averaging the (squared) lengths of all other links that do not touch the 
additional vertices. For all three considered measures this expectation value stays 
rather constant, (g^) ~ 5.1, independent of the cutoff/. This was also found for the 
completely regular triangulation. 

Fig. 1 shows (qi) as a function of / and one realizes that the heigth of the spikes 
increases with decreasing /. By means of numerical methods one cannot decide 
whether (qi) diverges in the limit / ^ or not. However, (qi) decreases with 
smaller a and the results show clearly that in the scale invariant limit, a ^ 0, the 
tendency to develop spikes at the inserted irregularities is suppressed. Thus, we 
use the 'nearly scale invariant' measure, a = j^^ also for computations with non- 
vanishing gravitational coupling and apply the lowest limit on the fatness, / = 10~^. 

Fig. 2a displays the average link length (qi) at the inserted vertices as a function 
of (3 in comparison to the expectation value (g^) of the regular part. At /3 > for 
every data point I.IM sweeps have been performed and after 100k thermalization 
steps every lOO*'' sweep has been recorded and used for the average values. A rather 
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interesting behavior is observed: As long as /3 < /3i ~ 0.025 a phase is found with 
monotonically increasing values of (qi) while (g^) stays almost constant (in fact it 
slightly decreases for (3 < 0.02). At /3i both expectation values suddenly jump to 
significantly larger values and enter between /3i < /3 < /32 ~ 0.04 a region with 
{qi) ~ 10"^ and (g^) ~ 10. At /32 the system performs a second jump into an ill- 
defined phase with diverging observables, known from the computations on regular 
hypercubic triangulations, which one can resolve for a finite cut off / > only. 

However, since the rotation to the Euclidean sector is not uniquely defined for 
quantum gravity, the sign of the coupling parameter (3 is not fixed a priori [3, 6]. 
We therefore present results also for negative gravitational couplings, /3 < 0, in 
Fig. 2a. Decreasing (3 stepwise, one observes a well-defined phase with homogeneous 
configurations, (qi) ~ (qr)- Thus, only 200k sweeps after 50k thermalization steps 
have been performed and measurement has been taken every lO*'' sweep for each data 
point. At /3o ~ —0.11 a transition to an ill-defined region with collapsing simplices 
occurs. Since the lattice size is small, it does not make sense to try to determine 
the order of this phase transition and we postpone this question to further studies. 

To obtain further insight we plot the expectation value ((^g) of the fatness as 
a function of (3. This observable (3) does not distinguish between 4-simplices con- 
taining irregular vertices. Nevertheless, all three transitions are refiected as seen in 
Fig. 2b. 

An important physical quantity is the expectation value (R) of the curvature 
as a function of (3. In Fig. 2c we observe an increase at (3i ~ 0.025 where the link 
lengths at the extra vertices become rather large. It is followed by another increase 
at (32 ~ 0.04 where the remaining regular vertices tend to develop spikes. This 
second transition represents the well-known transition from the well-defined to the 
ill-defined phase of the regularly triangulated Regge theory on hypercubes. Turning 
to negative gravitational couplings we find a decrease of (R) at /3o ~ —0.11 to an 
ill-defined phase with growth of the link lengths at both types of vertices. 



4 Higher-order terms for non-regular Regge lattices 

The introduction of an additional R^ term was initially proposed to cure the problem 
of the unbounded Regge action. However, there is numerical evidence that a phase 
with well-defined finite expectation values exists for simplicial quantum gravity even 
without such a term. Moreover, it has been claimed that an additional higher-order 
term ensures that the transition from the well-defined phase to the region of 'rough 
geometries' is of 2"'' order, accompanied by large correlation lengths, indicating a 
non-trivial continuum limit [7]. 

At this point it is important to remember that the Regge approach deals with a 
large number of unphysical degrees of freedom [5, 16]. On the triangulated hyper- 
cubic 4-torus one has 15 link lengths per vertex in contrast to the 10 independent 
components of the metric tensor. Notice that 4 of these 10 components can be fixed 
by choice of coordinates and the remaining 6 degrees of freedom are still subject of 
constraints refiecting the Bianchi identities. Respecting coordinate invariance, i.e. 
the diffeomorphism group, within the Regge calculus would lead to a non-trivial 
and non-local measure [15]. But until now, numerical studies within the Regge 
approach have ignored these difficulties and this might be justified far from the con- 
tinuum limit with the diffeomorphism symmetry completely broken. However, if one 
reaches the continuum limit, this question becomes important. (For a description 
of these problems in two dimensions see Refs. [16, 17].) Therefore, it is interesting 
to study higher-order gravity on general, non-regular triangulations with varying 
coordination numbers and degrees of freedom. H a continuum limit can be reached, 
recovering the diffeomorphism group, the universal behavior of the Regge approach 
should not depend on the choice of the incidence matrix. 

We thus now study the infiuence of the lattice construction in the presence of 
an additional term —aJ2t^t^t^t~ i^ ^^^ action with the 4- volume per triangle Vt 
given by barycentric subdivision [7]. Following the literature we set a = 0.005 where 
a 1** order transition at (AfSt) ~ was ruled out in favor of a 2 order transition 
[7]. To allow for a comparison with the results of hypercubic triangulation we use the 
uniform measure, a = X = 1. Our computation completely parallels those of Ref. [7] 
with the minor exception of posing a lower limit on the fatness, (fy^ y f = 10~^, and 
the major difference of using triangulations with Nq = 31,34 and 40 vertices, i.e. 



the regular and two non-regular triangulations with 3 and 9 additional vertices. 

In Fig. 3 we present results for the curvature on the regularly triangulated lattice, 
No = 31, which seems to behave as described in the literature [7]. At /3 ~ 0.25 the 
system reaches (R) ~ and a comparison of our data shows good agreement in 
the well-defined region, (R) < 0, /3 < 0.25. Different to the literature, we found 
stable expectation values even for (R) > and would locate the phase transition at 
(3 > 0.28. However, this is of minor significance since the lattice size is small, but 
indicates that it will be difficult to determine the order of the phase transition even 
for regular triangulations. 

Considering on the other hand the non-regular triangulations one observes that 
the additional action term suppresses the spikes at the inserted vertices. Fig. 4a 
compares directly the expectation value (qi) of the link lengths at the irregular ver- 
tices with {qj.) averaging over the links of the regular part for Nq = 34. Surprisingly, 
(qi) stays below (g^) for small /3, and we reduced in all i?^-term simulations the 
number of iterations per data point performing 50k sweeps to reach equilibrium and 
200k sweeps to gather statistics while measuring expectation values every lO*'' step. 
As in the case without R^ terms, first the extra vertices tend to build spikes at 
/3i ~ 0.16 followed by the remaining vertices at (32 ~ 0.25. 

The expectation value ((^g) averaged over the whole lattice does not distinguish 
between the different types of vertices, but refiects the behavior of (qi) and (g^)- 
Fig. 4b suggests a transition at /3i ~ 0.16 from the vertices with low coordination 
number and indicates a second transition to occur at (32 ~ 0.25 with all 4-simplices 
involved. 

The 3 additional vertices lead to a shift in the average curvature towards positive 
values and as seen in Fig. 4c (R) jumps at (3i ~ 0.16 from (R) ~ 2.6 to (R) ~ 40.3. 
A preliminary computation with 9 inserted vertices, Nq = 40, exhibits a further shift 
towards larger curvature. A comparison of Fig. 4c with 4a shows that the transition 
at (3i affects only the additional vertices whereas the transition at (32 from the regular 
vertices leads to an increase of (R) that is unlimited in principle for / ^ 0. 

To summarize, the effect of vertices inserted via b.s.d. steps is similar as in 
the pure Regge action discussed in the previous section. The additional higher- 
order term only shifts the critical values (3i,(32. Since the lattice size is small in our 



study, one has to be aware of the preliminary character of our findings. However, the 
results seem to be of importance for future simulations within the Regge approach for 
random triangulations on larger lattices. In general, for non-regular triangulations 
with varying local coordination numbers the infiuence of the incidence matrix on 
the phase transition seems to be essential: Even with an additional R^ term the 
possible existence of a 2 order phase transition at i? ~ 0, as observed for regular 
triangulations, is covered by another transition that occurs at vertices with low 
coordination number. 

5 Conclusion 

The Regge calculus provides a direct route to systematic approximations of the 
quantum-gravity path-integral. In a general setting one would investigate the path 
integral on random triangulations with arbitrary topologies. As a first step in this 
direction we considered non-regular triangulations by inserting additional vertices 
into a regularly triangulated 4-torus using barycentric subdivisions. 

We find that for vanishing gravitational coupling the inserted vertices develop 
spikes, as observed in earlier studies. Different to the regular, hypercubic triangula- 
tion, the infiuence of the measure becomes non-negligible, such that the formation 
of spikes is suppressed in the scale invariant limit. 

Using a nearly scale invariant measure we observe for negative gravitational cou- 
plings, > /3 > /3o ~ —0.1, a well-defined phase characterized by homogeneous 
configurations with the spikes at the irregularly inserted vertices strongly damped. 
If this property holds for arbitrary random triangulations one should investigate 
this phase more detailed, especially the character of the transition at /3o. For pos- 
itive couplings we observe two transitions, at /3i ~ 0.025 and /32 ~ 0.04, indicated 
by sudden changes of the average link lengths at the extra and regular vertices, 
respectively. The spikes at the vertices with low local coordination numbers are 
significantly larger than the average link length. 

The introduction of an additional higher-order R^ term in the action with 
a = 0.005 strongly suppresses these spikes. However, the critical behavior of the 
system seems to differ from that observed on regular, hypercubic triangulations. As 
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for pure gravity, the additional irregularities undergo a transition, at /3i ~ 0.16, in 
addition to that observed for the regularly triangulated part at (32 ~ 0.25. In future 
work, the phase structure of the Regge theory has to be examined for general trian- 
gulations and two main questions have to be dealt with: Do random lattices exhibit 
several transitions? Do the (two or more) transitions merge on larger lattices? 

Acknowledgment 

This work was supported in part by "Fonds zur Forderung der wissenschaftlichen 
Forschung" under Contract P9522-PHY. Part of the computations have been per- 
formed at the workstation cluster in SCRI at FSU, Tallahassee. W. Beirl thanks B. 
Berg and W. Bock for helpful and encouraging discussions. 

References 

[1] Ch. Misner, K. Thorne, and J. Wheeler, Gravitation (Freeman, San Francisco, 
1973). 

[2] T. Regge, Nuovo Cimento 19 (1961) 558. 

[3] J. Hartle, J. Math. Phys. 26 (1985) 804; J. Math. Phys. 27 (1986) 287; J. Math. 
Phys. 30 (1989) 452. 

[4] A. Nabutovsky and R. Ben-Av, Comm. Math. Phys. 157 (1993) 93; 
S. Catterall, J. Kogut, and R. Renken, hep-lat/9403019; 
J. Ambj0rn and J. Jurkiewicz, hep-lat/9405010; 
B. de Bakker and J. Smit, hep-lat/9405013. 

[5] M. Rocek and R. Williams, Phys. Lett. 104B (1981) 31; Z. Phys. C21 (1984) 
371. 

[6] B. Berg, Phys. Rev. Lett. 55 (1985) 904; Phys. Lett. B176 (1986) 39. 

[7] H. Hamber, in: Critical Phenomena, Random Systems, Gauge Theories, Pro- 
ceedings of the Les Houches Summer School, Session XLIII, Vol. 1, eds. K. 
Osterwalder and R. Stora (North Holland, Amsterdam, 1986) p. 375; Phys. 
Rev. D45 (1992) 507; Nucl. Phys. B400 (1993) 347. 

11 



[8] W. Beirl, E. Gerstenmayer, and H. Markum, Phys. Rev. Lett. 69 (1992) 713; 
W. Beirl, E. Gerstenmayer, H. Markum, and J. Riedler, Phys. Rev. D49 (1994) 
5231. 

[9] B. Berg, B. Krishnan, and M. Katoot, Nucl. Phys. B404 (1993) 359. 

[10] W. Beirl, E. Gerstenmayer, H. Markum, and J. Riedler, Nucl. Phys. B (Proc. 
Suppl.) 30 (1993) 764; 

W. Beirl, E. Gerstenmayer, and H. Markum, in: Quantum Physics and the 
Universe, Proceedings of the International Symposium, eds. M. Namiki, Y. 
Aizawa, K. Maeda, and I. Ohba (Pergamon, Oxford, 1993) p. 605; 
W. Beirl, Ph. D. thesis (TU-Wien, 1993), unpublished. 

[11] Kiihnel and Lassmann, Discrete Comp. Geom. 3 (1988) 169. 

[12] J. Cheeger, W. Miiller, and R. Schrader, Comm. Math. Phys. 92 (1984) 405; 
Led. Notes in Phys. 160 (Springer, Berlin, 1982) p. 176. 

[13] G. Feinberg, R. Friedberg, T. D. Lee, and H. Ren, Nucl. Phys. B245 (1984) 
343. 

[14] M. Bander, Phys. Rev. Lett. 57 (1986) 1825. 

[15] A. Jevicki and M. Ninomiya, Phys. Rev. D33 (1986) 1634. 

[16] W. Bock and J. Vink, hep-lat/9406018. 

[17] J. Ambj0rn, B. Durhuus, and J. Fr6hlich, Nucl. Phys. B257 (1985) 433. 



12 



100 



10 - 

















1 




- 














n 


(7=1 '- 


'. Q- ■ 


■ ■□. . 












A 


(7 = 1 : 


■ 




■Q-. 


■ ■□. 











- = i^ - 


A. 


■A. 






■q- . 


■ ■□• . 












■■A. 


■■A. 


■■A. 


■A.. 


■A.. 


■■A. 




■ 0- 
.... 1 


■ o- ■ 


. .0. . 


■ ■0- ■ 


■ -lOi 




■ -0 ■ 


■ 0- ■ 
1 


■G....O : 






. . . . 1 




■ -0 ■ 



10 



100 



1000 



/[lo- 



Fig. 1. Dependence of the simulation on the measure for vanishing gravitational coupling. A 
non-regular triangulation of the 4-torus based on the construction of Kiihnel and Lassmann 
[11] with 3 additional vertices obtained by bary centric subdivision has been considered. For 
each value of the measure parameter a = 1,1, and jq the lower limit of the fatness is 
decreased stepwise by a factor -^ towards 10~®. The average next-neighbor distance of the 
additional vertices diverges for / ^ as long as a is significantly larger than 0. Notice the 
logarithmic scale. Statistical error bars are in the size of the symbols. 
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length (qi) at the additional vertices and (qr) of the regularly triangulated part are shown 
in plot (a) in a logarithmic scale. The system undergoes transitions in the positive coupling 
region at /3i ~ 0.025 where the extra vertices develop spikes and at /32 ~ 0.04 where the rest 
follows. For negative couplings the well-defined phase continues with homogeneous config- 
urations, (qi) ~ (ir), with the spikes at the additional vertices suppressed. At /3o ~ —0.11 
a transition to an ill-defined region with collapsing simplices is observed. The expectation 
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as a function of /3. The system reaches (R) ~ at /32 ~ 0.25. Our data is compared to 
R((]) = -A{(]c - PY with A = 3.794, /3e = 0.2443, 6 = 0.624 as reported in [7]. 
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Fig. 4. Simulation with additional i?^ term. The average squared link length [qi) at the 
additional vertices and [q^) of the regularly triangulated part are shown in plot (a) in a 
logarithmic scale. At (3i ~ 0.16 a transition of the link lengths at the additional vertices is 
observed and at (32 ~ 0.25 the remaining links start to grow. The expectation values of the 
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